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Abstract. We propose a model for the acceleration of charged particles in interplanetary space that appear 
during quiet time periods, that is, not associated with solar activity events like intense flares or coronal mass 
ejections. The interaction of charged particles with modeled turbulent electromagnetic fields, which mimic the 
fields observed in the interplanetary medium, is studied. The turbulence is modeled by means of a dynamical 
system, the Gledzer-Ohkitani-Yamada (GOY) shell model, which describes the gross features of the Navier- 
Stokes equations. The GOY model is used to build a 3-D velocity field, which in turn is used to numerically 
solve the ideal magneto-hydrodynamic (MHD) induction equation, while the electric field is calculated from the 
ideal Ohm's law. Particle acceleration in such an environment is investigated by test particle simulations, and 
the resulting energy distributions are discussed and compared to observations of suprathermal electrons and ions 
during quiet periods in interplanetary space. 
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The observations of high energy, suprathermal particle 
populations in interplanetary space reveal the presence 
of a rich var i ety o f physical characteristics and processes 
I see IKeame 2 GBs) for a review!. The information recov- 
ered from time profiles of particle fluxes, energy spec- 
tra, element abundances, ionization states, etc., is essen- 
tial in determining many properties of the sources and 
of the mechanisms of particle acceleration. The remark- 
able heterogeneity found in the detected energetic par- 
ticle events is related, in large part, to the existence of 
different sources of acceleration in interplanetary space. 
Suprathermal populations can be produced by solar flares, 
collisionless shock waves driven by Coronal Mass Ejections 
(CMEs), Corotating Interaction Regions (CIRs), or the 
heliospheric termination shock. 

For many years, the attention of researchers has been 
focused on suprathermal particles with energies above a 
few hundred keV, due to the fact that the spacecraft in- 
struments lacked the sensitivity needed to investigate the 
energy range from the solar wind thermal plasma up to 
a few hundred keV. In recent years, however, with the 
launch of the Ulysses, Advanced Composition Explorer 
(ACE), and WIND spacecrafts, this gap has been filled 
and observations of "quiet time" suprathermal particles. 
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that is, not associated with the abrupt energization events 
mentioned above, have become available. 

The electron spectrum, measured between ~5 eV and 
~100 keV by the 3-D Pla sma and Energetic Particles 
Experiment ( Lin et alJl995t) on the WIND sp acecraft dur - 
ing a quiet period, has been investigated by 'LirJ l)l998ll . 
A Maxwellian core dominates the spectrum from ~5 eV 
to ~50 eV, while a hotter popu lation, the so-called solar 
wind halo i Feldman et al.lll975(l . takes over in the range 
between ~100 eV to ~1 keV, due to the escape of coronal 
thermal electrons with temperature of 10^ K. However, 
these WIND observations have made possible the iden- 
tification of a third, much harder component, which has 
been denoted the "super-halo" , with energies from ~2 keV 
up to ^100 keV and an approximate power law shape 
with exponent ~ 2.5. The angular distribution of these 
"super-ha lo" electrons is nearly isotropic. According to 
iLirJ l|l998l) . this high energy population is not solar in ori- 
gin, since this would imply a continuous production and 
escape of electrons with such energies from the Sun. It 
has been suggested ()Linlll998l) that the "super-halo" tail 
is due to acceleration by CIRs beyond 1 AU, but clear 
evidence for correlations with CIRs or solar active regions 
have not yet been found. 

The velocity distributions of solar wind ions from 
0.6 to 100 keV/e, measured using the Solar Wind 
Ion Composition Spectrometer (SWICS) instruments on 
Ulysses and ACE (,Cloeckler et al...l992..1995.b have been 
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Studied in iGloeckleil l)l999|) . iGloeckler et~al] l)2000|) . and 
iGloeckleil (1200,11) . One of the most important findings of 
these works is that the speed distributions of H+ , He"*" and 
He ions show well developed, approximate power law 
tails during quiet time periods, that is, far from shocks, 
CIR compressions and other disturbances. This indicates 
the presence of a population of highly suprathermal ions at 
all times. The power law exponents of these tails, which 
extend over the whole measurement range, are between 
~ 5 and ~ 5.5 in the slow, in-ecliptic solar wind, and 8 
in the super-quiet fast wind coming from polar coronal 
holes. Comparing ACE observations at ^ 1 AU to Ulysses 
observations at ~ 5 AU, the authors also found that the 
tails are continuously regenerated in the out-flowing so- 
lar wind, overcoming the cooling related to the wind ex- 
pansion. The main question arising from these observa- 
tions is how these ubiquitous suprathermal ions are pro- 
duced in the quiet solar wind when there are no shocks. 
CIRs o r other disturbances observed locallv. lTe Roux et al.l 
suggested that pickup ions might be accelerated 
by large-scale turbulent electric fields directed along the 
background magnetic field. In order to explore this pos- 
sibility, they presented a numerical model for gyrotropic, 
pitch-angle dependent pickup ion transport between the 
Sun and the Earth based on standard kinetic theory for 
charged particles. The ion kinetic transport equation used 
in their model includes a Gaussian random value repre- 
sentation of the large-scale field-aligned electric field fluc- 
tuations averaged over the characteristic length and time 
scales of MHD turbulence in the low-latitude solar wind. 
The authors choose the standard deviation of these elec- 
tric fields by requiring the reproduction of observed ac- 
celerated pickup ion spectra and in this way they were 
able to obtain a qualitative agreement of the results of 
their model with the suprathermal He+ spectra in the 
slo w low-latitude solar wind, observed at ~ 1 AU. 

iKirsch fc Mall Q003) recently presented an analysis of 
interplanetary suprathermal ions, based o n measurements 
perfo rmed with the SMS experiment I Gloeckler et alJ 
Il995i) on the WIND spacecraft, in the range 0.5-225 keV/e. 
The authors investigated particle bursts in which high en- 
ergy protons in the range 5-100 keV are observed in asso- 
ciation with distinct decreases of the magnetic field mag- 
nitude. They considered only events not associated with 
shocks, CIRs, and magnetospheric disturbances. As a con- 
sequence of the the observed magnetic field behaviour, 
they suggest that these bursts could be the result of a 
local reconnection process, or, alternatively, they propose 
that inductive electric fields [i.e. V x £^ = —{l/c)dB/di\ 
could be a possible explanation for the observed particle 
acceleration. 

Here, we investigate with a new approach the pos- 
sibility that the high energy particles observed in the 
interplanetary space during quiet time periods are due 
to a process of stochastic acceleration in the turbulent 
electromagnetic fields present in the heliospheric plasma. 
Recently, the stochastic acceleration process in turbulent 
electromagnetic fields has been investigated with numeri- 



cal experiments in which test particle simulations are per- 
formed in field configurations that are obtained from the 
solution of the magneto-hydrody namic (MHD) equations 
l|Nodes et alJ2003HDmitruk et aL2003^) . The authors sug- 
gest that these simulations could potentially be applied to 
astrophysical problems. As an alternative approach, other 
authors have studied the motion of test particles in elec- 
tromagnetic flelds built up by means of suitable models 
for particular app lications, like the Ea rth's magnetotail 
dVell^ri et alJll99'8l) or the solar corona ijArzner & Vlaho^ 
2004) . With respect to other approaches, test particle sim- 
ulations offer the possibility to describe some peculiar 
features of particle acceleration in turbulent fields, espe- 
cially the possibility for particles to be trapped, possibly 
in or around strong, coherent electric field regions, lead- 
ing to effective acceleration even for low init ial energies 
l|Ambrosiano et al.lll98llDmitruk et al.ll2003|) . 

In the present paper, we present a model for 
stochastic, quiet time particle acceleration in the in- 
terplanetary space, based on turbulent electromagnetic 
fields constructed by means of the so- called G l edzer- 
Ohkitani-Yamada (GOY) shell model llGledzeil Il97.'^ 
lYamada k, Ohkitanilll98H) . which mimics the nonlinear 
dynamics of fiuid turbulence. The use of such a simplified 
description of turbulence implies that we can describe only 
some basic features of the nonlinear interactions occurring 
in turbulent fluids, with the advantage that we can simu- 
late turbulent electromagnetic fleld conflgurations without 
the computational difficulty of solving the full MHD equa- 
tions. The magnetic field is determined through the MHD 
induction equation, assuming weak magnetic fields so that 
the back reaction onto the plasma can be neglected. The 
electric field is given by Ohm's law. In this framework, we 
perform test particle simulations, concentrating mainly on 
the energetics of the injected particles, i.e. electrons and 
ions. 

A basic description of the model is given in Sect. [3 In 
Sect 121 we present the numerical simulations performed 
and the results obtained from them. Discussions and con- 
clusions are given in Sect. 0] 



2. The model 

The model describes the acceleration of charged test 
particles in turbulent electromagnetic fields, as obtained 
from a dynamical system model of turbulence. In gen- 
eral, the macroscopic description of a plasma is given 
by the ma gneto-hydrodynamic (MHD) equations [see e.g. 
iBovd fc"Sanderson (2003) ]. Here, we consider particle ac- 
celeration events occurring in regions of the interplanetary 
space where the magnetic field is weak, so that we can 
restrict ourselves to weakly magnetized plasmas, which 
implies that the temporal evolution of the plasma is gov- 
erned by the velocity field. In order to build up the ve- 
locity field configurations, we use the so-called COY shell 
model, through which it is possible to generate a turbu- 
lent, incompressible 3-D velocity field V(r,i), as will be 
explained in Sect. 12.11 
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The velocity field is used to obtain numerical solutions 
of the MHD induction equation for a perfectly conducting 
plasma, namely 



dB 

'dt 



Vx{V xB) 



(1) 



where B is the magnetic field. The dissipative term /iV^B 
(where is the magnetic difFusivity) is not taken into ac- 
count, since the plasma in the interplanetary space can be 
considered coUisionless to a good approximation [see e.g. 
[Montgomery ( 1983 )]. The electric field is then computed 
from the ideal Ohm's law 



E = --V X B 



(2) 



where the resistive term rjj (with rj being the resistivity 
and j the current density, respectively) is again neglected. 
As already mentioned before, the feedback of the mag- 
netic field on the velocity field (that is, the effect of the 
Laplace force) is not considered, as the velocity field is 
given independently and the MHD momentum equation 
is not considered. In other words, we suppose that the 
magnetic energy density is much smaller than the kinetic 
energy density of the turbulent flow. 

To investigate the acceleration of test particles in the 
electromagnetic fields generated as described above, we 
consider the relativistic equations of motions of a charged 
particle in an external electromagnetic field: 



dp 

lit 



dr 

dt 



= qE +-V X B, 



(3) 



(4) 



where r, v, and p are the position, velocity, and momen- 
tum of the particle, respectively, c the speed of light, and 
q the charge of the particle. For the numerical integration, 
we express Eq. Q in terms of velocity (777,71; = p) , 



^ = ^E- 
dt -fm 



jnic 



-V X B 



"fmc 



Mv ■ E), (5) 



where 7 = 1/ -y/l — w^/c^, and m is the mass of the parti- 
cle. 

2.1. Construction of the 3-D fluid velocity field 

The 3-D, time-dependent velocity field V{r, t) is con- 
structed by means of the so -called GOY shel l model for 
turbulence. Shell models [see iBohr et all l|2nn(ll) for a com- 
plete review] are dynamical systems designed to repre- 
sent in a simplified way the spectral form of the equations 
which de scribe turbulent fluids. They were origina l ly pro - 
posed by'Obukhov' (" iQTl^ . lDesnvanskv fc Novikovl l)l974|) . 
and iGle dzcr ( IQT^i) for hydrodyn a mic turbulence. The 
GOY shcU model llGledzerl Il973t lYamada fc Ohkitanil 
has been extensively investigated, both analitically 
and numerically (,Yamada fc Ohkitanii.1987; .Jensen et alJ 



Il99lt iBiferale et alJll995|) . In the following, we describe 
some basic characteristics of the GOY model. 

The main idea of the GOY shell model is to mimic the 
Navier-Stokes equations by a dynamical system in which 
the velocity fleld fluctuations at different length scales are 
represented by scalar variables u„(i). To this aim, the 
Fourier space is divided into N shells, with the associated 
wave number denoted by A;„, where the shell index n is 
discrete. The scalar, complex variable Un{t) is associated 
with the 77-th shell, and the nonlinear dynamics of tur- 
bulent fluids are modeled by quadratic couplings among 
nearest and next nearest neighbour shells, following the 
assumption that the nonlinear interactions are local in the 
k space. The coefficients of the nonlinear terms are deter- 
mined by imposing the conservation of the ideal invariants 
of the Navier-Stokes equations. The e quations of evolut ion 
of the dynamical variables u„(i) are (jBohr et al.ll200Cl() 

d ,2 

dt^^'- 



7fc„_iU*_iM*_2) + Smnfn , 



where n • 



(6) 



The parameter v stands for the kinematic viscosity, while 
Smnfn IS a stochastic forcing term acting on the shell m, 
one of the first shells, providing a constant average energy 
fiux into the system {6mn is the Kronecker symbol). The 
wave numbers are chosen to follow the relation 



koh" 



(7) 



where fco and h are constant (/i > 1, usually h = 2), and 
n > 1. The shells are thus equally spaced in a logarithmic 
scale, which is justified by the fact that in fully developed 
turbulence the energy spectrum in the nonlinear, inertial 
range follows a power law. 

One of the main advantages of shell models over nu- 
merical simulations of the Navier-Stokes equations is that 
they can be investigated at much higher Reynolds num- 
bers. They provide a good description of the scaling prop- 
erties of fully developed turbulence in the inertial range, 
even if, being scalar models, they do not include informa- 
tion about the spatial structures of turbulence. From the 
scalar variables u„(t) we can generate an incompressible 
velocity field V{r,t) in the r eal 3-D space by applying 
a simple numerical algorithm l|Bohr et al.ll2000|) . The dy- 
namical variables w„ {t) of the shell model are supposed to 
represent the coefficients of a Fourier expansion with wave 
vectors fc in a shell of radius |fc| — kn. We introduce a set 
of vectors fc„: 

fen — knS-n , (8) 

where Sn = {e^ \ ei^\ el^''} are randomly chosen vectors 
of unit norm. The components Vj{x, t), j = 1, 2, 3, of the 
velocity field are obtained by the analogue to an inverse 
Fourier transform. 



N 

y,(r,<) = VcWK(i)e^'^"-^ 



n=l 



+ C.C. 



(9) 
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where the coefficients C„ are of order 0(1). 

In order to satisfy the incompressibility constraint 
V ■ V = 0, the vectors — {en \ en \ ei^''} and the coef- 
ficients Cn'' must satisfy the condition 

3 

^CWe(f)=0, Vn. (10) 

3. Numerical simulations and results 

The numerical investigation consists of two main steps: 
1) the calculation of the magnetic and electric fields from 
Eqs. and (0), in which the velocity field given by Eq. @ 
is used; 2) test particle simulations through the numerical 
solutions of Eqs. © and ((SJ). 

3.1. Magnetic and electric field calculation 

The GOY shell model equations Eq. 10 are solved by us- 
ing a 4-th order Runge-Kutta integration algorithm, using 
N = 22 shells. The kinematic viscosity in the shell model is 
assumed tohev = 10^^. Once the shell model has reached 
a statistically stationary state, we start the numerical inte- 
gration of the MHD inducti on equation [Eg. dH)], by using 
the Wilson upwind scheme l)Hawlev et al.ll9 84). imposing 
free outflow boundary conditions and using the velocity 
field V{r, t) from Eq. I©. The initial magnetic field Bo(r) 
is given by a random perturbation constructed through a 
sum of Fourier modes Bo{r) = J2k -^o(fc) cos(fc -r + ipf^), 
with Gaussian distributed amplitudes |Bo(fc)|, random 
phases (pf^ and with the constraint k Bo{k) = imposed, 
as it follows from V • B = 0. The size of the grid is 64'^. 
The evolution of the system is followed over a time inter- 
val 2re, where is a typical eddy turnover time, given 
by Te = L/V„-ns , L being the size of the simulation do- 
main and VJ-ins the rms velocity. The condition that the 
kinetic energy density is much smaller than the magnetic 
energy density is checked during the time evolution. We 
also monitor the value of V • B, evaluated through a stan- 
dard finite differences scheme, and verify that it does not 
vary significantly with respect to the initial value, which 
can be considered zero within the numerical error. 

The equations are solved in non-dimensional form, and 
suitable rescaling factors are applied to describe quiet 
time periods in interplanetary space. The applied forc- 
ing term leads to an rms velocity field intensity of Kms — 
3.2 X 10^ cm s~^. This value is slightly larger that the 
rms velocity estimates obtained from a 30 year dataset 
of solar wind o bservations, which gave values around 
7 X 10^ cm s^i ifBreech et al.ll2003 '). This means that in 
our simulations we are assuming a slightly enhanced tur- 
bulence level with respect to the average. The linear size of 
the simulation box is assumed to be L = 2.2 x 10^" cm. The 
rms value of the magnetic field intensity, after the rescal- 
ing, is Sims — 7.4 X 10~^ G, while the rescaled rms value 
of the electric field intensity is i?rms — 5.7 x 10~^ stat- 
volt cm~^. These values are slightly larger than the rms 



values obtained from long time datasets of solar wind ob- 
servations. 

To illustrate the 3-D structure of the magnetic and 
electric fields, in Fig. we present a 3-D visualization of 
the regions where the magnetic field and electric field in- 
tensities exceed the thresholds Bth = 1.2 x 10^'' G and 
Eth = 1.2 X 10~^ statvolt cm~^, respectively. 

In Fig.|2we present the Probability Density Functions 
(PDFs) of the three electric field components, collected at 
a fixed time from the entire simulation box. In order to 
compare the different PDFs, the variables are first trans- 
lated to zero mean and then normalized to their standard 
deviation, so that all the PDFs have zero mean and unit 
standard deviation. As it can be seen, these PDFs are not 
Gaussian, they exhibit clear exponential tails. This result 
is in qualitative agreement with the one-point PDFs of 
the observed interplanetary indu ced electric fields (lEF) 
presented inlBreech et all l|2003|) . In a more recent work 
teor riso-Valvo et al.'2004\. it has been shown that the sta- 
tistical properties of the IFF depend on the wind velocity, 
an effect that we do not model here (see discussion in Sect. 
EI). 

3.2. Test particle simulations 

In this subsection, we present some results obtained from 
test particle simulations of electrons and ions in the elec- 
tromagnetic field configurations generated as described 
in the previous subsection. Specifically, the magnetic and 
electric field configurations obtained at the time 2Te dur- 
ing the evolution of Eq. are used. The particle motion 
equations ^ and ^ are solved with a 4th order Runge- 
Kutta, adaptive step-size scheme. The magnetic and elec- 
tric field configurations are kept constant during the time 
wc monitor the particles, assuming a much slower evolu- 
tion time for the fluid velocity fleld than for the test par- 
ticles. Since the magnetic and the electric fields are given 
only at a discrete set of points, both fields are interpolated 
with a local, 3-D linear interpolation to provide the field 
values in between grid-points, wherever they are needed 
for the integration scheme. The initial time step used for 
the integration is set to O.ltg, where tg is the gyration 
period of the particles at the starting point. 

The charged test particles are injected at random posi- 
tions within the simulation box, with velocities extracted 
from the tail of a Maxwellian distribution with tempera- 
ture T = 1.15 X 10^ K (corresponding to ^ 10 eV), which 
represents a typical value in the interplanetary space at 
1 AU from the Sun. The threshold velocity used to select 
only the particles in the tail of the Maxwellian is 2vth, 
where vth is the thermal velocity. The choice to use only 
the tail of the initial Maxwellian distribution is due to the 
fact that wc do not model the coUisional processes occur- 
ring in the interplanetary space plasma, so we assume that 
only the high energy part of the distribution participates 
in the acceleration process. If a particle leaves the simu- 
lation box before the end of its tracing time interval, it is 
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Fig. 1. 3-D visualization of the magnetic field intensity | 
of the electric field intensity \E\ above the threshold Eth 
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Fig. 2. PDFs of the three components of the electric field: 
the solid, dashed, and dot-dashed lines represent the PDFs 
of Ex , Ey , and Ez respectively. 



reinjected into the box at a random point on the surface 
opposite to the one through which it had left. The maxi- 
mum tracing time interval used in this work is t = 300 s, 
which is much smaller than the typical collision times for 
electrons and ions with velocities larger th an twice the 
thermal velocity [see e.g. iMontgomervTljlQSSl) ] . so that we 
can neglect coUisional effects. 

In Fig. 121 a sample trajectory of a typical test electron 
is shown. It can be seen that, during its irregular motion, 
the particle visits both regions where it remains trapped 
for some time and regions where its motion exhibits long 
"jumps" , as could be expected in a turbulent field envi- 
ronment. In Fig.^ we report part of the time evolution of 



I above the threshold Bth = 1.2x10 "'G (left panel) and 
: 1.2 X 10"^ statvolt cm (right panel). 




Fig. 3. Trajectory of a test electron. 



the kinetic energy of a typical test electron (upper panel) , 
the associated electric field intensity along the trajectory 
(middle panel) , and the cosine of the angle a between the 
particle velocity and the electric field (lower panel). 

In Fig. 13 a sample trajectory of a test proton is shown. 
As is expected, the typical shape of the proton trajectories 
is substantially different from the one of the electrons, due 
to the much larger mass of the protons. In Fig.|B| we report 
part of the time evolution of the kinetic energy of a test 
proton (upper panel), the associated electric field intensity 
along the trajectory (middle panel), and the cosine of the 
angle a between the particle velocity and the electric field 
(lower panel). From this figure, it is more clear, with re- 
spect to the figure referring to electrons (Fig. ^ , that the 
largest kinetic energy variations are associated both with 
intense electric field spikes and with extended regions of 
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Fig. 4. Part of the time evolution of the kinetic energy 
of a test electron (upper panel), the associated electric 
field intensity along the trajectory (middle panel), and the 
cosine of the angle a between the particle velocity and the 
electric field (lower panel). 



almost coalignment between particle velocity and electric 
field. 

In Fig. [71 we report the kinetic energy PDFs of 10'* test 
electrons at the initial time i = s, and at three successive 
times t = 1 s, t = 30 s, and t = 300 s respectively. At 
t = Is, the evolution of the initial Maxwellian distribution 
produces a tail, extending to ~ 5 keV, which does not show 
a clear, unique power law form, but is more of a double 
power-law shape. As the time increases, the high energy 
tail tends to become exponential with maximum energy 
up to ~ 100 keV 

In Fig. |S1 we report the kinetic energy PDFs of 10^ 
test protons at the same times as for the electrons. The 
evolution of the proton distributions is qualitatively sim- 
ilar to the electrons', although the PDF tail at t = 1 s 
now displays a reasonably clear power law scaling that 
persists until 30 s. However, for t — 300 s the tail is nar- 
row and steep, reminiscent of both a steep power-law and 
an exponential distribution, with maximum energy up to 
~ 100 keV. 




Fig. 5. Trajectory of a test proton. 

Since observations are also available for other ions, and 
in particular for He ions, we also investigated the kinetic 
energy PDF of He+ ions. The time evolution of this PDF, 
shown in Fig. is quite similar to the proton case. For 
t = 1 s, the tail of the PDF extends to lower energies 
than for the protons, again being a power-law shape that 
persists until i = 30 s. For t — 300 s, energies up to 
^ 100 keV are reached, with a narrow and steep tail that 
is difficult to classify. 

4. Discussion and conclusion 

In this paper, we have presented a model for turbulent 
particle acceleration based on a dynamical system descrip- 
tion of turbulence. The aim of the model is to describe 
long lasting particle acceleration processes occurring in 
the turbulent interplanetary space during quiet time peri- 
ods, which lead to the appearance of suprathermal tails at 
all times in the energy distributions of electrons and ions 
(extending up to ^ 100 keV). 

The acceleration process has been investigated by per- 
forming test particle numerical simulations in the electro- 
magnetic fields obtained by numerically solving the ideal 
MHD induction equation, which is driven by the velocity 
field that is calculated using a dynamical system model 
(the so-called GOY shell model) of turbulence. This ap- 
proach implies that the magnetic fields are supposed to be 
weak, or that the magnetic energy density is much smaller 
than the kinetic energy density of the flow. 

The PDFs of the electric field components have been 
shown to be non-Gaussian, exhibiting exponential tails. 
The presence of roughly exponential tails in the one-point 
PDFs of the interplanetary induced electric fields (lEF) 
has been recently shown based on data analysis performed 
on 30 years of me asurements that w ere acquired by dif- 
ferent spacecra fts (Breech et al."2003). However, more re- 
cently l|Sorriso-Valvo et al. ,2004') . it has been shown that 
using homogeneous datasets with respect to wind veloc- 
ity and solar activity, that is, considering short datasets 
and separating the data according to slow and fast wind 
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Fig. 6. Part of the time evolution of the kinetic energy of 

a test proton (upper panel), the assoeiated electric field 
intensity along the trajectory (middle panel), and the co- 
sine of the angle a. between the particle velocity and the 
electric field (lower panel). 



streams, the exponential tails are recovered only in the 
radial components of the electric field, that is, the compo- 
nent along the Sun-Earth direction, which coincides with 
the mean wind direction. Our model is obviously not able 
to reproduce in detail these statistical properties of the ob- 
served lEF, since we do not take into account some basic 
features of the solar wind structure, e.g. its mean bulk ve- 
locity, mean magnetic field structure, etc. In other words, 
we consider only the effect of field fluctuations related to 
the presence of turbulence. 

With our approach, we have been able to obtain basic 
physical insights into the process of particle acceleration 
due to turbulent electric fields in weakly magnetized plas- 
mas, and to investigate the possibility that this mecha- 
nism plays a role in the acceleration of charged particles 
to suprathermal energies in interplanetary space during 
quiet periods. The trajectories obtained from the simula- 
tions indicate that the particles alternately visit regions 
in which they are trapped for some time, due to the ef- 
fects of turbulence, and other regions where long "jumps" 
are observed. The observed motion of the test particles 
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suggests that a Fokker-Planck description of the diffusion 
process underlying the particle acceleration is not suffi.- 
cient to achieve a complete characterization of the prob- 
lem. This is one of the main reasons why a test particle 
approach is adequate in studying such situations. 

We found that, starting from an initial thermal pop- 
ulation, both in the case of electrons and ions, the initial 
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Fig. 9. Kinetic energy probability density function of 10^ 
test He+ ions at the initial time t = (dotted curve), and 
after t — \ s (dashed curve), t = 30 s (dot-dashed curve), 
t = 300 s (solid curve). 

Maxwellian energy distributions evolve in time, giving rise 
to power-law tails for shorter times, which become very 
steep and narrow for the longest time we monitor the 
particles, so that it is difficult to discriminate between 
exponential and power-law distributions. At the maxi- 
mum time we allow, the particles reach energies up to ^ 
100 keV. The fact that our model is able to reproduce the 
observed energies suggests that a stochastic acceleration 
mechanism resulting from the turbulence that is devel- 
oped in the interplanetary space can be at the origin of 
the ubiquitous suprathermal tails observed during quiet 
time periods, even if the detailed shape of the observed 
distributions, which exhibit approximate power law tails, 
is not very well recovered by our model. 

The difficulty in reproducing the observations in detail 
at this stage is related to some limitations of our model in 
its current form. Our main simplifying assumptions are: 
(1) We have assumed that the magnetic energy density is 
much smaller than the kinetic energy density of the fluid. 
This assumption is often not fulfilled in the interplanetary 
space, due to the presence of a strong Alfvenic compo- 
nent i n the solar wind turb ulence [see e.g. lColdstein et iiD 
||1995(1 : ITu fc MarschI l(l995^) ]. (2) The detailed p roperties 
of intermitt e ncv in solar wind turbulence fsee e.g.lBurlagal 
()l99lLll992|) : lMa,rsch fc iTil l)l 993^ : Fa,rbone et a,1.l (|l flfl.'j) ] 
are not included in the model, and this could have an ef- 
fect especially on the high energy part of the kinetic energy 
PDFs. (3) The resistivity, which modifies the small-scale 
structure of the electric field, is treated as a constant, not 
taking its possible d y namic evolution i i ito ac count [see 
Inmitnik et ^ hm± : lArzner fc V1a,hoj l)2004ll ]. (4) The 
large scale magnetic field structure, which was not in- 



cluded in the present study, but which nevertheless would 
be interesting to be investigated, will have a less deciding 
influence on the energetics. 

In order to overcome the limitation (1), a dynamical 
system for MHD turbulence must be considered, instead 
of a hydrodynamical model, whereas for the limitation (2), 
a more appropriate description of intermittency should be 
introduced in the reconstruction of the 3-D fields from the 
1-D scalar fields that are yielded by the shell models. The 
presented results suggest that on modifying our model in 
this way, better compatibility with the observations can 
be reached. 
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